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Shock  wave  interactions  in  hypervelocity  flow 

S.  R.  Sanderson  and  B.  Sturtcvant 

California  Institute  of  Technology,  Pasadena,  CA  91125,  USA 

Abstract.  The  impingement  of  shock  waves  on  blunt  bodies  in  steady  supersonic  flow  is  known 
to  cause  extremely  high  local  heat  transfer  rates  and  surface  pressures.  Although  these  problems 
have  been  studied  in  cold  hypersonic  flow,  the  effects  of  dissociative  relaxation  processes  are  un¬ 
known.  In  this  paper  we  report  a  model  aimed  at  determining  the  boundaries  of  the  possible 
interaction  regimes  for  an  ideal  dissociating  gas.  Local  analysis  about  shock  wave  intersection 
points  in  the  pressure-flow  deflection  angle  plane  with  continuation  of  singular  solutions  is  the 
fundamental  tool  employed.  Further,  we  discuss  an  experimental  investigation  of  the  nominally 
two-dimensional  mean  flow  that  results  from  the  impingement  of  an  oblique  shock  wave  on  the 
leading  edge  of  a  cylinder.  The  effects  of  variations  in  shock  impingement  geometry  were  visualized 
using  differential  interferometry.  Generally,  real  gas  effects  are  seen  to  increase  the  range  of  shock 
impingement  points  for  which  enhanced  heating  occurs.  They  also  reduce  the  type  IV  (Edney  1968 
a,b)  interaction  supersonic  jet  width  and  influence  the  type  II— III  transition  process. 

Key  words:  Shock-on-shock  interaction.  Shock  impingement,  Hypervelocity  flow,  Dissociation, 
Relaxation,  Heat  transfer 

1.  Introduction 

Shock  impingement  phenomena,  which  inhibit  the  further  development  of  hypervelocity  vehicles, 
are  discussed  in  the  literature  by  many  authors,  notably  by  Edney  (1968  a,b).  For  the  inviscid, 
compressible  flow  of  a  perfect  gas  a  sufficient  set  o  of  dimensionless  parameters  to  describe  any 
quantity  in  the  flow  is,  <p  =  <p[M,  7,  fi\ ,  A,  I'],  where  M  is  the  free  stream  Mach  number,  7  is  the  ratio 
of  specific  heats,  /3,  is  the  impinging  shock  wave  angle,  A  describes  the  position  of  the  impinging 
shock  wave  with  respect  to  the  body  and  T  is  a  set  of  parameters  defining  the  body  geometry. 
Consider  the  case  of  a  given  gas  and  fixed  free  stream  condition,  impinging  shock  strength  and 
body  geometry.  The  only  remaining  dependence  is  then  the  location  of  the  impingement  point 
relative  to  the  body,  i.e.,  <*>  =  #[A].  As  the  incident  shock  wave  is  translated  relative  to  the  body 
it  potentially  intersects  with  al!  possible  shock  strengths  of  both  positive  and  negative  slope.  On 
the  basis  of  an  experiment  such  as  this,  with  spherical  and  modified  spherical  bodies,  Edney  (1968 
a,b)  observed  and  categorized  six  interaction  regimes;  known  as  types  I-VI.  Edney  rationalized  the 
observed  flow  fields  through  local  analysis  about  shock  wave  intersection  points  in  the  pressure- 
flow  deflection  angle  plane  (p-6  plane).  The  key  conclusion  drawn  from  such  an  analysis  is  the 
role  of  the  three-shock  solutions  (or  A-shocks)  and  this  is  discussed  further  in  §2.  Based  on  the 
assumption  of  straight  shocks  in  the  vicinity  of  the  interaction,  conclusions  were  drawn  about  the 
global  nature  of  the  flow.  When  the  impinging  shock  wave  is  weak,  the  flow  downstream  of  a 
A-shock  is  characterised  by  low  Mach  number  flow  behind  the  strong  transmitted  wave  and  high 
Mach  number  flow  downstream  of  the  weak  reflected  shock.  On  dimensional  grounds,  the  inertia  of 
the  supersonic  stream  is  dominant.  Thus,  the  flow  field  for  the  global  type  IV  interaction,  whereby 
a  supersonic  jet  penetrates  a  region  of  low  subsonic  flow  [see  Edney  (1968  a,b)  or  Fig.  4],  can  be 
solved  approximately  as  a  free  streamline  flow  up  to  some  unknown  length  scale.  Typically  this 
length  scale  is  specified  in  terms  of  the  width  of  the  jet. 

According  to  Edney ’s  model,  the  heat  transfer  is  determined  by  the  attachment  at  the  body 
surface  of  the  shear  layers  generated  at  the  shock  impingement  points.  A  correlation  was  obtained 
between  local  pressure  and  local  heat  transfer  rate  at  the  surface.  A  further  observation  made 
by  Edney  is  that  the  jet  curvature  increases  as  the  intersection  point  moves  up.  Some  recent 
contributions  to  the  literature  are  noted  in  the  references  of  this  paper. 

2.  Local  analysis  for  an  ideal  dissociating  gas 

The  flow  in  smooth  regions  and  across  shock  discontinuities  is,  in  principle,  understood.  The 
predominant  concern  is  then  points  where  surfaces  of  discontinuity  intersect.  Local  analysis  about 
intersection  points  in  the  p-b  plane  provides  a  useful  tool  for  analysis.  This  technique  is  well 
known  in  the  Mach  reflection  literature — see  for  example  Hornung  (1986).  Our  purpose  here  is  to 
generalize  these  tools  to  the  rase  of  oblique  translational  shock  waves  where  both  the  upstream  and 


downstream  conditions  are  non-equilibrium  states  with  respect  to  the  internal  degrees  of  freedom 
of  the  gas.  This  analysis  is  performed  for  a  one-dimensional  shock  and  subsequently  extended  to 
oblique  waves. 

Consider  the  dimensionless  form  of  the  continuity  and  momentum  equations  that  apply  across 
a  one-dimensional  translational  discontinuity  and  throughout  the  downstream  relaxation  zone, 

1  p-  1 

pu=l,  p=\  +  (1,2) 

where  P  =  p\/p\U\*.  The  symbols  have  the  usual  meanings  and  generally  <p  = 

We  use  Lighthill’s  (1957)  model  for  the  thermodynamics  of  an  ideal  dissociating  gas  (IDG): 

P  =  2~(1  +«)p7' ,  &  =  2~[(4  +  a)7'  +  a0J  ,  (3,4) 

where  k  is  Boltzmann’s  constant,  m  is  the  mass  of  one  atom  of  the  gas,  Bd  is  a  temperature 
characterizing  the  dissociation  energy  and  a  is  the  dissociated  mass  fraction.  Using  equation  (4) 
the  dimensionless  energy  conservation  equation  becomes, 


rn  r  • 

(4  +  0'i)~  +  <*2  +  rj  =  Ho  , 
D  p2 


(5) 


where  the  dimensionless  parameters  are  D  —  0d/T\ ,  l\  =  m>j\/kdd  and  H0  =  2mho/k9d.  H0  is 
the  conserved  stagnation  enthalpy. 

Importantly,  for  a  nonequilibrium  binary  mixture,  three  parameters  are  sufficient  to  define  the 
thermodynamic  state  of  the  upstream  gas.  Here  we  specify  P,  H-d  and  aq.  The  remaining  param¬ 
eters  are  then  given  by  the  following  identities,  obtained  from  H 0  written  in  terms  of  conditions 
upstream  of  the  discontinuity, 


Hp  -  ay  _  4  -f  Q) 

1  +  2^1^  ~  H0-K-ai  ' 


(6,7) 


From  equation  (3)  it  follows  that  T  -  and,  eliminating  T  and  D  from  equation  (5),  with 

the  pressure  given  by  (2),  leads  to  the  following  quadratic  equation  for  the  density, 


(//0-o2)p2- 


^+A>/>  +  i^  =  o, 

'  p 


A 


=  2  KP 


4  +  02 
1  +  02 


(8) 


One  solution  describes  the  variation  of  density  downstream  of  a  translational  shock.  The  remaining 
solution  describes  relaxation  of  the  (possibly)  nonequilibrium  upstream  state  without  any  discon¬ 
tinuity.  Both  solutions  are  parameterized  in  terms  of  the  dissociation  mass  fraction  o2,  whose  rate 
of  reaction  need  not  yet  be  specified.  This  interpretation  for  the  case  O]  =  0  is  due  to  Hornung 
(1979). 

Equation  (8)  may  be  combined  with  an  equilibrium  expression  to  obtain  the  asymptotic  state 
far  downstream  from  the  translational  shock.  For  the  equilibrium  IDG, 
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P±e-D/t 
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(9) 


This  leads  to  a  second  transcendental  equation  relating  «2  and  p.  Numerically  solving  equation 
(8)  with  o2  =  given  by  (9)  gives  the  equilibrium  shock  solution  for  specified  upstream  con¬ 
ditions  P,  Ho  and  o  t.  Further,  solving  equation  (8)  with  n->  =  o;  gives  the  frozen  solution  valid 
immediately  downstream  of  the  translational  discontinuity. 

The  extension  to  oblique  shock  waves  is  obtained  in  the  usual  way.  The  normal  components, 
of  l>  and  I\  are  given  in  terms  of  ff  the  angle  of  the  oblique  shock  wave  by 


l< i ».  =  A’i  siir  .1 


(10,11) 
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These  are  used  in  (6),  (8),  (9)  and  (2),  and  the  the  flow  deflection  angle  6  is  determined  from 
p  =  tan  /?/  tan(/7  -  6). 

To  facilitate  the  computation  of  multiple  shock  jumps  for  the  interaction  problem  we  must 
find  a  convenient  means  of  determining  P  and  Ho  downstream  of  the  oblique  shock  Trivially 
H  o3  =  Ho,.  It  can  be  sh  that  A'2  =  /\'i  (cos2  /?  +  sin2  0/p2).  Pi  then  follows  from  (6). 

3.  Continuation  of  s  r  solutions 

While  the  p-6  plane  r<  ,  .nation  is  useful  for  interpretation  of  solutions  it  is  difficult  to  deal 
with  the  case  distinctions  that  arise  as  the  number  of  parameters  is  increased.  One  technique  for 
reducing  the  amount  of  information  is  to  consider  only  the  singular  cases  that  delineate  regions  of 
solutions  with  similar  forms.  Various  special  cases  arise  that  are  significant  in  terms  of  transition 
to  Mach  reflection.  Discussion  of  the  significance  of  these  points  can  be  found  in  the  literature 
[e.g.,  Hornung  (1986)].  Now  generalize  these  cases  to  the  asymmetrical  interaction  of  oblique  shock 
waves; 

i.  Interaction  of  shock  waves  that  produce  a  Mach  stem  with  zero  curvature.  This  is  a  general¬ 
ization  of  the  von  Neumann  condition  in  Mach  reflection. 

ii.  Sonic  flow  downstream  of  the  transmitted  shock  waves. 

iii.  Coalescence  of  the  strong  and  weak  regular  solutions  such  that  the  shock  loci  intersect  only 
tangentially.  This  is  a  generalization  of  the  maximum  deflection  condition  in  Mach  reflection. 

iv.  Sonic  flow  downstream  of  the  incident  shock  waves. 

v.  One  transmitted  shock  wave — A-shocks. 

Given  the  solution  for  one  such  set  of  these  singular  points,  solutions  for  neighboring  values  of  the 
parameters  may  be  determined  using  path  following  techniques  [see  Keller  (1987)]  for  either  frozen 
or  equilibrium  chemistry.  The  frozen  solution  applies  in  regions  near  the  shock  wave  intersection 
point  and  the  equilibrium  solution  is  valid  far  from  the  shock  waves.  The  regular  4-shock  intersec¬ 
tions  map  into  the  2nd  and  4th  quadrants  of  Fig.  1  and  lie  on  the  surface  bounded  by  the  axes  and 
the  von  Neumann  curves.  Outside  this  region  the  only  solutions  are  the  A-shocks  which  exist  only 
along  lines.  If  the  impingement  point  on  a  blunt  body  maps  to  a  region  of  Fig  1  where  no  regular 
solution  exists  then  the  flow  must  deform  in  such  a  way  that  the  A-shock  condition  is  maintained 
at  the  intersection  point.  This  highlights  the  physical  mechanism  determining  the  gross  nature  of 
the  flow  field. 

It  is  interesting  to  consider  various  classical  problems  in  this  framework,  a)  Shock  detachment 
from  a  wedge;  <r i  =  (To.  Hornung  and  Smith  (1979)  observed  significant  non-equilibrium  effects. 
This  is  consistent  with  the  separation  of  the  frozen  and  equilibrium  sonic  lines  in  Fig.  1.  b)  Mach 
reflection;  cj\  —  —a2.  Hornung  et  af.  (1979)  observed  that  non-equilibrium  effects  were  slight. 
This  is  consistent  with  the  proximity  of  the  frozen  and  equilibrium  von  Neumann  lines  in  Fig.  1 . 
c)  Shock  impingement  problem;  <rt  =  const.  Little  experimental  evidence  exists  for  high  enthalpy 
effects  in  this  case.  For  weak  impinging  shock  waves,  shown  by  the  vertical  line  in  Fig.  1,  the 
separation  of  the  frozen  and  equilibrium  A-shock  solutions  is  relatively  large.  This  suggests  some 
effect  on  the  transition  from  type  II  to  type  111.  This  occurs  in  the  vicinity  of  the  sonic  line  which 
is  itself  sensitive  to  non-equilibrium  processes. 


4.  Experimental  results 

We  have  conducted  an  experimental  investigation  of  the  nominally  two-dimensional  mean  flow 
that  results  from  the  impingement  of  an  oblique  shock  wave  on  the  leading  edge  of  a  cylinder. 


This  problem  is  abstracted  from  the  complex  configurations  which  occur  in  the  flow  about  flight  ->r _ 

vehicles  and  contains  much  of  the  physics.  These  initial  experiments,  comprising  35  shots,  were  0^ 

conducted  in  the  GALC1T  T5  free  piston  reflected  shock  tunnel  with  nitrogen  test  gas  at  nozzle 
reservoir  enthalpies  of  3  MJ/kg  and  12.5  MJ/kg  and  reservoir  pressures  of  12  MPa  and  25  MPa, 

respectively.  The  reservoir  gas  was  expanded  through  a  contoured  axisymmetric  nozzle  with  area  rj _ 

ratio  109.5  to  yield  a  nominal  Mach  number  of  5.5  in  the  test  section.  The  cylinder  was  37.5  mm  _ - 

in  diameter  with  aspect  ratio  4.5  and  the  flow  was  deflected  7  degrees  by  the  incident  planar  shock 
wave,  'fhe  effects  of  variations  in  shock  impingement  geometry  were  visualized  using  infinite  fringe  . 

differential  interferometry  with  vertical  beam  shear.  A  Q-switched,  frequency  doubled,  NdYAG  . - 

laser  was  used  as  the  light  source.  _  ty  OffdaS 
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sgn(/?2)  (n/2-\p2\)/(n/2-fjJ) 


a,  =  sgn(j?,)  (n/ 2-|/?,|)/  n/2-^J 

Fig.  1.  Map  of  interaction  regimes  for  two  intersecting  oblique  shock  waves  of  angles  3[  and  /?2- 
Shock  angles  are  normalized  so  that  (T  =  0  corresponds  to  a  normal  shock  wave.  Mach  waves  correspond 
to  a  =  ±1.  The  free  stream  conditions  are  taken  as  H o  =  0.5,  P  =  0.03  (or  =  5  with  */  =  4/3), 
Otl  =  0  and  =  10s. 


Figs.  2-4  illustrate  type  II,  III  and  IV  interactions  at  stagnation  enthalpies  of  3MJ/kg  and 
12.5  MJ/kg.  The  incident  shock  is  concave  down  at  the  edge  of  test  section  and  some  fringe  shift 
is  observed  below  this  wave.  This  complicates  interpretation  of  the  images.  However,  the  main 
centerline  disturbance  is  generally  at  the  top  of  the  line-of-sight-i nt egrat ed  image. 

At  the  lowest  impingement  point  reported  in  this  paper  the  type  I— 1 1  transition  has  already 
occurred  with  the  incident  shock  wave  impinging  somewhat  below  the  lower  sonic  line  (Fig.  2).  In 
the  low  enthalpy  case  a  Mach  stem  connects  the  A-shock  pattern  at  the  impingement  point  to  the 
bow  shock  of  the  cylinder.  At  the  intersection  point  a  strong  vortex  sheet  is  generated  separating 
the  tipper  region  of  subsonic  flow  from  the  lower  region  of  superset  ic  flow.  No  transmitted  wave  is 
resolvable  at  the  upper  end  of  the  Mach  stem  which  appears  to  join  smoothly  into  the  bow  shock 
in  the  vicinity  of  the  sonic  line.  The  upper  portion  of  the  bow  shock  remains  symmetrical.  In  the 
high  enthalpy  case  the  bow  shock  remains  undisturbed.  However,  a  point  of  inflection  is  observed 
immediately  above  the  intersection  point. 
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(h)  12.5  MJ/kg 


interactions. 


In  the  typo  III  experiments  (Fig.  3)  the  incident  shock  wave  intersects  the  how  shock  wave  in 
t  lie  subsonic  region,  somewhat  below  the  geometric  stagnation  point.  For  both  enthalpies  the  Mach 
stem  has  merged  with  the  bow  shock  wave  and  the  A-shock  pattern  is  preserved  at  the  intersection 
point.  'The  influence  of  the  impinging  shock  wave  is  now  global  and  both  the  radius  of  curvature 
and  the  standoff  distance  of  the  asymmetrical  upper  portion  of  the  bow  shock  has  increased  to 
match  the  A-shock  pattern  at  the  impingement  point.  Dominated  by  the  inertia  of  the  supersonic 
stream,  the  shear  layer  is  deflected  only  in  the  immediate  vicinity  of  the  cylinder  (c.f.  discussion 
in  Ijl ).  Generally,  tie  standoff  distance  decreases  with  increasing  enthalpy  and  this  influences  the 
impingement  of  the  shear  layer  on  the  body.  In  the  high  enthalpy  case,  the  resulting  compression 
waves,  that  are  obs<  rved  as  a  series  of  nearly  vertical  white  lines  below  the  impingement  point, 
steepen  the  lower  shock.  Kinks  are  observed  where  these  waves  meet  the  supersonic  continuation 
of  the  bow  shriek.  \\  ak  reflected  expansions  are  observed  as  horizontal  disturbances  emanating 
from  the  intersection  points.  Note  the  analogy  between  this  phenomenon  and  the  mechanism  of 
complex  Mach  reflection. 

When  the  inri.l-  nt  shock  wave  impinges  in  tie'  vicinity  of  the  geometrical  stagnation  point 
a  type  IV  flow  result'  (Fig.  -I  i  file  A-shock  pattern  is  preserved  and  the  supersonic  transmitted 
portion  of  the  bow  shock  forms  a  second  inverted  A-shock  with  the  subsonic  continuation  of  the 
bow  shock  below  tie  interaction,  l  he  flow  behind  the  oblique  shock  wave-  connecting  the  two  A- 
sliocks  is  supersonic  and  forms  a  jet  that  turns  upwards  under  the  action  of  the  pressure  gradient 
produced  by  the  sec.  aid  A-shock  'I  he  strongly  curved  shear  lay.  r  generated  at  the  upper  A-shock 
forms  large  pi  nines  as  it  passes  above  the  cylinder.  I  lie  standoff  distance  and  jet  width  are  reduced 
in  the  high  enthalpy  case  but  this  also  depends  strongly  on  tie  location  of  the  impinging  shock. 
\  portion  of  the  sup-  fsonic  f)..w  is  now  turned  so  that  it  pass,  s  above  tie  cylinder  and  hence  the 
stagnation  sireamliu-  must  pass  through  the  supersonic  jet.  I  lie  stagnation  density  is  higher  for 
streamlines  that  pass  'liroiigli  t  le- supersonic  jet  than  for  st  ream  lines  that  cross  the  adjacent  strong 


(a)  3  MJ/kg  (h)  1*2.5  MJ/kg 


Fig.  4  l>iir<-r<'iil ial  iiiterferograms  of  1  y p«-  IV  interactions. 

bow  shock.  The  higher  density  and  large  unsteady  velocity  gradients  produced  by  the  impingement 
of  the  supersonic  jet  provide  the  mechanism  for  local  increase  of  the  heat  transfer  rate 

5.  Conclusions 

bocal  analysis  about  shock  wave  intersection  points  in  the  //-/>  plane  has  been  extended  to  the  ideal 
dissociating  gas  model.  The  application  of  continuation  methods  to  singular  solutions  illustrated 
the  role  of  A-shooks  and  determined  the  inflm  nee  of  real  gas  effects  on  local  analysis.  Kxperiments 
conducted  in  a  free  piston  shock  tunnel  provided  insight  into  real  gas  effects  on  the  global  flow. 
Generally,  real  gas  effects  are  seen  to  increase'  the  range  of  shock  impingement  points  for  which 
enhanced  heating  occurs.  They  further  influence  the  Mach  stem  form  in  type  II  Hows,  the  deflection 
of  the  shear  layer  in  type  III  flows  and  the  length  scales  associated  with  the  type  IV  jet. 
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